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Co NMR evidence for charge ordering below Tco ~ 51 i^ in Nao.5Co02 
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The C0O2 layers in sodium-cobaltates Naa;Co02 may be viewed as a spin S = ^ triangular-lattice 
doped with charge carriers. The underlying physics of the cobaltates is very similar to that of the 
high Tc cuprates. We will present unequivocal ^ Co NMR evidence that below Tco ~ 51 A", the 
insulating ground state of the itinerant antiferromagnet Nao.5Co02 (Tjv ~ 86 K) is induced by 
charge ordering. 

PACS numbers: 71.27.+a, 71.30.+h, 76.60.-k 



The discovery of unconventional superconductivity in 
sodium-cobaltate Nai/3Co02[H20]4/3 (T^ ~ 4.5 K) 
[Hi 0,9 has generated major excitement in the condensed 
matter community. Co ions in the C0O2 layers take 
a mixed valence state 00+"*^^, and form a triangular- 
lattice. Since €0"*+ and 00"^+ ions nominally have spins 
of S" = 2 and S = 0, respectively, one may consider 
the C0O2 layers as a S" = ^ triangular-lattice doped 
with charge carriers, in analogy with the doped Cu02 
square-lattice in the high Tc cuprate superconductors. 
The fingerprints of Co spins are everywhere in NaxCo02: 
Nao.82Co02 is an itinerant antiferromagnet (Tjv ~ 21 K) 
[J|; Nao.7Co02 is a "Curie Weiss metal" with large para- 
magnetic susceptibility @; Nan.5Co02 is an itinerant an- 
tiferromagnet (Tjv ~ 86 iiT) 0, and undergoes a myste- 
rious metal-insulator transition below ~ 51 ivT [^. 

In the case of the cuprates, doped holes tend to segre- 
gate themselves from the underlying S = ^ square-lattice 
to form static charge stripes [6| . Whether the charge or- 
dering in stripes is related to or competing against the 
mechanism of high Tc superconductivity is controversial 
[7|. A major question we address here is whether or not 
the cobaltates also exhibit a similar phenomenon in the 
triangular-lattice geometry. The insulating ground state 
of Nao.5Co02 below ~ 51 if has attracted considerable 
attention because of speculation that it may be charge 
ordered [1, B H, HI. However, earher ^SNa NMR mea- 
surements in Nao.5Co02 showed no evidence of the emer- 
gence of a charge ordered state below ^ bl K |lll.[l3.[l3j. 

In this Letter, we take a more direct approach to prob- 
ing the charge environment in C0O2 layers, by measuring 
the EFG {Electric Field Gradient) of Co sites with zero- 
field 5^ Co NMR. The EFG tensor is the second derivative 
of the Coulomb potential, and hence is directly related 
to the local charge density. We will present unequivocal 
evidence that the insulating ground state of Nao.5Co02 
below Tco ~ 51 iiT is indeed the consequence of charge 
ordering. 

A major distinction between cobaltates and cuprates 
is that the Na+ ions in cobaltates can spatially or- 
der for certain values of Na concentration, x. Be- 
sides donating electrons, ordered Na"*" ions exert a pe- 



riodic Coulornb potential on the C0O2 triangular-lattice 
[i[il[il,[ii,[il[3. in the case of Nao.5Co02, electron 
and neutron diffraction measurements have suggested 
that Na+ ions form zigzag chains [a, [iJ, ^M,- This re- 
sults in the presence of two structurally inequivalent Co 
sites, each forming a chain structure[l5| as shown in Fig. 
1. One of the nearest neighbor sites of Co(l) is occu- 
pied by a Na+ ion, while Co(2) has no Na"*" ions in the 
nearest neighbor sites [1^. Since the Coulomb poten- 
tial from Na+ chains attracts electrons, the valence of 
Co(l) sites, Co"'''^'^"'', is smaller than that of Co(2) sites, 
(2;q+3.5+(5 rp-j^g average valence of Co ions in Nao.5Co02 
is +3.5. The valence of Co(2) sites is closer to Co+^ 
[S ~ i), hence Co(2) exhibits strong spin fluctuations 
above T/v ~ 86 K [T^, and develops a sizable ordered 
magnetic moment ~ 0.26/iB within the C0O2 plane be- 
low Tat [3, [201 ■ On the other hand, Co(l) sites have 
much weaker spin fluctuations above T/v [3 ■ The upper 
bound of the ordered moment at Co(l) sites is as small 
as ~ 0.04^B at 8 K and too small to be detected by po- 
larized neutron scattering techniques [^. ^^Co zero field 
NMR is a very powerful technique for detecting small hy- 
perfine magnetic field Bhf, which is proportional to the 
ordered magnetic moments. Yokoi et al. took advantage 
of this high sensitivity, and demonstrated the presence of 
small ordered moments at Co(l) sites along the crystal 
c-axis below T/v ^ i& K . This led them to propose two 
possible spin structures as shown in Fig.l 20]. 

Although the existence of multiple Co sites in 
Naa;Co02 is sometimes referred to as a consequence of 
charge ordering persisting up to room temperature, it 
is important to realize that these distinct behaviors of 
Co(l) and Co(2) sites in Nao.5Co02 are directly linked 
with the periodic Coulomb potential arising from Na^ 
zigzag chains. Moreover, charges on the Co sites are mo- 
bile, and Nao.5Co02 is metallic above Tco ~ 51 K [^. 
In contrast, strong electron-electron correlation effects in 
high Tc cuprates induce a self-organized pattern of doped 
carriers in the form of charge stripes 6] . The spatial peri- 
odicity of the charge ordered state is different from that 
of the Coulomb potential of the underlying lattice. It 
is this type of self-organized charge pattern that we are 



searching for in Nao.5Co02. 

In Fig. 2, we present typical zero- field ^^Co (nuclear 
spin I ~ ■^) NMR lineshapes from a piece of electro- 
chemically deintercalated Nao.5Co02 single crystal (mass 
^ 25mg) I2ll. 122] ■ The observed lineshapes are similar 
to Yokoi'sjil. Below Tat - 86 K, we observe 7 NMR 
peaks of Co(2) sites for transitions between nuclear spin 
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and /» 



2m- 



■i, where integer —3 < m < +'3. 



In the frequency range above 5 MHz, we have also suc- 
cessfully detected Co(l) zero- field NMR signals for the 
m = to m = -1-3 transitions. The intensity of NMR 
signals decreases in proportion to the square of the fre- 
quency, and measurements below 5 MHz are formidable. 
We summarize the temperature dependence of the ob- 
served peak frequencies fm of Co(l) sites in Fig. 3. A 
remarkable feature of Fig. 2 and Fig. 3 is that all 4 peaks 
of Co(l) sites split into two precisely below Tco ^ 51 K. 
In contrast, Co(2) sites show no anomalies. In what fol- 
lows, we will refer to the two inequivalent Co(l) sites 
as Co(la) and Co(lb). In passing, the comparison be- 
tween the zero-field and low-field NMR lineshapes indi- 
cates that the NMR signals near 6 MHz, shown in green, 
are not related to Co(l) sites (1^. We tentatively at- 
tribute these unidentified signals to defect sites caused 
by slight deviation of Na concentration from 0.5 and/or 
minor disorder. We also confirmed that both Co(la) and 
Co(lb) lines split in weak magnetic field applied along 
the c-axis J13j. This means that both Co(la) and Co(lb) 
sites have up-spins as well as down-spins. 

To understand the effects of the EFG on NMR line- 
shapes, we need to theoretically compute the resonance 
frequency fm of the m-th transition as fm = (-E(2m+i)/2 — 
E(2m-i)/2)/h, where the energy levels En are eigenvalues 
of the standard ^^Co nuclear spin Hamiltonian, 
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The first term of the Hamiltonian represents the nuclear 
quadrupole interaction [2^ between the EFG and the 
^^Co nuclear quadrupole moment. The diagonalized ten- 



sor of the nuclear quadrupole interaction {uq 



-'o) 



is proportional to the EFG tensor (dcf)'^ / dX"^ , d(jP /dY"^, 
dcf)'^ / dZ"^) , where cj) is the total Coulomb potential seen 
by ^^Co nuclei, and X, Y and Z represent three orthogo- 
nal, principal axes (by following the convention, we define 
From the comparison of zero-field 






< Wq\ < \vf 



NQR (Nuclear Quadrupole Resonance) results at 110 K 
(see Fig. 2) and high-field NMR, we found that the main- 
principal axis Z coincides with the crystal c-axis within 
experimental uncertainties (^ 5 degrees). This means 
that the other principal axes, X and Y , lie within the 
C0O2 plane. Poisson's relation for the Coulomb inter- 
action sets a constraint, Vq + i^q + i^q =0. Hence we 
have only two independent parameters in the first term 
of eq.(l), i'q (we now denote Z = c) and the anisotropy 
parameter 77 = {vq—Vq)/vq. ry is a measure of the devia- 
tion from axial symmetry of the EFG tensor with respect 



to the crystal c-axis. In the paramagnetic state at 110 K 
(> Ttv), we deduced from the ^^Co NQR lineshape that 
j/^ = 2.799(5) MHz and 77 = 0.300(1) for Co(l) sites, and 
i/^ = 4.040(5) MHz and r? = 0.400(1) for Co(2) sites. 

The second term in the Hamiltonian represents the 
Zeeman interaction between Co nuclear spins and the 
local magnetic field, B, at the position of the observed 
nuclear spin. 7„ = 27r x 10.054 MHz/Tesla is the ^^Co 
nuclear gyromagnetic ratio. In our zero-field NMR mea- 
surements, we apply no external magnetic field. Hence 
B = B/iy, where '^hf is the hyperfine magnetic field from 
ordered Co moments. B^j is proportional to the sublat- 
tice magnetization, M(r), of Co sites [20j. The ra = 
transition frequency depends primarily on the Zeeman 
term, /o - 7„|B/,/|. At 60 K, /o = 6.659 MHz and 19.780 
MHz for Co(l) and Co(2) sites, respectively, because the 
hyperfine magnetic field is |B;i/| ~ 0.6 Tesla at Co(l) 
sites, and ~ 1.9 Tesla at Co(2) sites. The separation be- 
tween adjacent peaks, {fm+i — fm), depends primarily 
on the EFG. For Co(l) sites, {fm+i - fm) -^ Vq -^ 2.8 
MHz, since the ordered moments point along the c-axis. 



For Co(2) sites, the separation is {fm+i — fm) 






'Q 



/2 ^ 2 MHz or less, because the ordered magnetic 
moments point within the ab-plane. The observed peaks 
fm are not evenly spaced, because the nuclear quadrupole 
and Zeeman interactions are comparable in the present 
case. This means that we must diagonalize the Hamil- 
tonian in eq.(l) exactly to fit the lineshapes. By relying 
on inaccurate second-order perturbation analysis, Yokoi 
et al. not only misidentified the peaks of Co(2) sites but 
were also unable to deduce the critically important infor- 
mation associated with the Co(l) sites |2(j |. 

Based on the exact diagonalization, we matched all the 
observed zero-field NMR peaks with theoretically calcu- 
lated values of fm as shown by vertical lines in Fig. 2. 
This allowed us to deduce the relevant NMR parameters 
of eq.(l), i.e. j/q, 77, the magnitude of -B^/, and the rel- 
ative polar angle {9, (j)) between Bhf and the c-axis. We 
found that r] shows no temperature dependence across 
51 K, and {9, cj)) change less than 2 degrees from the val- 
ues at 60 K (61 = 9 ± 1°, = 0°). Hence we focus our 
attention on the temperature dependence of Vq and Bhf 
of Co(la) and Co(lb) sites. 

Co(la) and Co(lb) have different values of i/q, as 
shown in Fig4b. This is direct proof that local charge en- 
vironment is different between Co(la) and Co(lb) sites. 
In principle, a differentiation of the EFG between Co(la) 
and Co(lb) sites could arise if a structural phase tran- 
sition doubles the unit cell of Na+ zigzag chains. To 
rule out such a scenario, we deduced Vq^No) at two 
structurally inequivalent Na sites, Na(l) and Na(2) [15j . 
from the measurements of -|-i to — ^ central transition 
and ±1 to ib^ satellite transitions of ^'^Na NMR in the 
same crystal |13| . As summarized in Fig. 4a, we found 
no hint of additional splitting or extra broadening for 
h'Q{Na) below Tco- The upper bound of the potential 
splitting of VQ{Na) at 4 K is less than ~ 1%, while 



the splitting of Vq between Co(la) and Co(lb) sites, 

Avq = Vq °^ '^Q " ' reaches as much as ~ 6%, as 
shown in Fig.4d. Therefore we conclude that Co(l) sites 
undergo a charge ordering transition at Tco ^51 K. 

The observed differentiation of Bhf between Co(la) 
and Co(lb) sites in Fig. 4c is also understandable as a 
consequence of charge ordering. Since the valence state 
is different below Tco ~ 51 K^ the number of electrons 
filling the t2g orbitals of Co(la) and Co(lb) sites will also 
be different. This naturally results in a different magni- 
tude of ordered moments at Co(la) and Co(lb) sites, 
and hence separate values of B^j. In Fig.4d, we plot 
the splitting of Bhf between Co(la) and Co(lb) sites, 

^Bhf = Bf^f "'' — Bf^f '. ABfif shows an identical 
temperature dependence as Ai^g within experimental un- 
certainties. This is not surprising, because to a good 
approximation, both A:^g and AB/jj should be propor- 
tional to the difference of the number of electrons occu- 
pying the Co 3d orbitals. Thus we may consider Avq 
and ABhf as the order parameter of charge ordering be- 
low 51 K. By fitting the temperature dependence down 
to 40 K {— O.STco) to a typical power law behavior, 
Aiy^ ~ ABhf ^ {T - TcoY' ■, we found the critical ex- 
ponent 13' = 0.3 ±0.1. Interestingly, this value is compa- 
rable to the critical exponent /3 = 0.28 ± 0.02 observed 
for the sub-lattice magnetization, M{T) ~ (T — Tm)^ at 
Tjv = 86 K for the antiferromagnetic Neel transition of 
Co(2) sites [9]. 

To conclude, we have presented unequivocal zero-field 
NMR evidence that static differentiation of charge den- 
sities develops below Tco ~ 51 K in the one dimensional 
chain structure of Co(l) sites. In contrast, Co(2) sites 



have only spin ordering and show no anomalies either 
in NMR lineshapes or vq across 51 K. In Fig.l, we 
show the possible charge ordering patterns with the 
smallest possible unit cell. Our observation of a striped 
configuration of the charge ordered state, consisting 
of Co(l) chains with charge and spin ordering and 
Co(2) chains with spin ordering only, is consistent with 
the two-fold symmetry observed for in-plane magne- 
toresistance by Balicas et al. Q. The splitting of Vq 
(and hence EFG) reaches as much as ~ 6 % between 
Co(la) and Co(lb) sites. We recafi that Vq at ^^Cu 
sites increases by ^ 10% in high Tc cuprates when the 
hole concentration x changes by 0.15 in the Cu02 plane 
of La2_2;Sr2;Cu04 [2j|. If we assume that the charge 



doping effects on the EFG are comparable between 
cobaltates and cuprates, we can crudely estimate the 
differentiation of the valence between Co(la) and Co(lb) 
sites as '^ 0.1. It is worth noting that t-U-V model 
calculations [lC| showed that a differentiation of valence 
of Co(2) sites as little as ~ 0.03 may be sufficient to 
drive C0O2 planes into insulating. On the other hand, 
none of the charge ordered patterns proposed within the 
existing theoretical frameworks [T^, [2g, [2^ are consistent 
with our experimental finding of charge ordering on 
Co(l) sites. Clearly, more theoretical works are required. 
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Fig.l. Ordered magnetic moments on Co(l) sites 
point either up (•) or down (x), while moments on Co(2) 
sites point within the plane (arrows) [20|. Co(l) sites 
have either in-plane antiferromagnetic order (left panel) 
or in-plane ferromagnetic order (right panel). In the 
latter case, Co(l) sites in the adjacent C0O2 layers have 
opposite spin orientation. Dark blue and light blue on 
Co(la) and Co(lb) sites, respectively, represent possible 
charge ordering patterns to be determined in this work. 
Red dashed lines represent the unit cell in each possible 
configuration. 



Fig.2. 59Co NMR fines at Co(l) sites (blue) and Co(2) 
sites (orange) at 110 K (paramagnetic state, Bhf = 0), 
60 K (Neel ordered state), 30 K and 4.2 K (charge 
ordered state). Vertical lines represent theoretical fit 
of resonance frequencies /,„ for Co(l) (black), Co(la) 
(dark blue), Co(lb) (light blue), and Co(2) (orange). 
Peak(s) near 6 MHz (green) are not associated with 
Co(l) or Co(2) sites (see main text). 

Fig.3. /„ (m = 0,1,2 and 3) for Co(l) (black) 
above Tco, and Co(la) (dark blue) and Co(lb) (light 
blue) below Tco- Vertical arrows roughly represent the 
magnitude of i^q for Co(la) and Co(lb) sites. 

Fig. 4. (a) Temperature dependence of VQ{Na) at 
Na(l) and Na(2) sites. The distribution of VQ{Na) 
is about the size of the symbols. (b) Temperature 
dependence of i^g, and (c) hyperfine field B^f in Co(l) 
sites. Color convention is the same as in Fig.2 and 
Fig.3. (d) /S.Uq : the splitting of Vq between Co(lb) and 
Co(la) sites deduced from Fig. 4b (left axis, diamond). 
^Bfif : the splitting of Bhf between Co(la) and Co(lb) 
deduced from Fig. 4c (right axis, filled circles). Solid 



curve shows a best fit to Av, 
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AB 
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{T - Tcof 

with a critical exponent /3' = 0.3. The dashed line shows 
Tco -^ 51 K. 
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